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ABSTRACT 

Recent observations have revealed the existence of enormously energetic ~ 10 61 erg 
AGN outbursts in three relatively distant galaxy clusters. These outbursts have pro- 
duced bubbles in the intra-cluster medium, apparently supported by pressure from 
relativistic particles and/or magnetic fields. Here we argue that if ^GcV particles 
are responsible then these particles are very likely protons and nuclei, rather than 
electrons, and that the 7-ray emission from these objects, arising from the interac- 
tions of these hadrons in the intra-cluster medium, may be marginally detectable with 
instruments such as GLAST and HESS. 
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1 INTRODUCTION 

Galaxy clusters are the 
bound systems in t he 
l|Giovannini fc Ferettil |2000| 



universe 
r 



largest 
In 



gravitationally 



add ition, radio 
Feretti et all 120041) and 
hard X-ray (|Rephaeli fc Gruberll2002l : iFusco-Femiano et al.l 
2004)) observations have revealed that a significant compo- 
nent of non-thermal particles can be found in such systems. 
The remaining tracer of non-thermal particles is high 
energy 7-ray emission, but no such signal has been firm ly 
detected from galaxy clusters so far l|Reimer et al.ll2003l ). 

Despite this non-detection a number of arguments sug- 
gest that galaxy clusters are potentially powerful emitters 
of hig h ener gy radiation. Volk. Aharonian fc Breitschwerdtl 
(1996) and iBerezinskv. Blasi fc Ptuskinl l|l997l ) recognised 
that hadronic cosmic rays (CRs) with energies of less 
than 10 15 eV accumulate within the cluster volume for 
the entire Hubble time. This CR component, together 
with the presence of target material in the form of 
the hot intra-cluster medium (ICM), will lead to very 
high energy (VHE) 7-ray production via i nelastic proton- 
proton collisions and subsequent n° dec ay (|Denniso n] |l980l ; 
IVolk, Aharonian fc Breitschwerdtl Il996l ). Furthermore, lep- 
tonic CRs are also capable of generating high energy elec- 
tromagnetic radiation: TeV electrons may up-scatter cos- 
mic microwave background (CMB) p hotons to 7-ray ener- 
gies in the inverse Compton processes (|Atovan fc VolklboOOl ; 
iGabici fc Bl"aslll2003l . liobl ). The lifetime of VHE electrons, 
however, is limited by IC and synchrotron losses to ~ 10 6 
years (for typical cluster magnetic field strengths). There- 
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fore, only recently injected electrons will contribute to the 
production of VHE 7-rays. Finally, if populations of ultra 
high energy (> 10 18 eV) CR protons exist in galaxy clus- 
ters, they will interact with CMB photons and produce elec- 
tron - positron pairs which will in turn radiate TeV photons 
(with a char acteristically hard energy spectrum) via the IC 
mechanism (jlnoue. Aharonian fc Su givama 2005). 

Several sources of CRs are plausible in galaxy 
clusters. Large scale shock waves caused by hierar- 
chical structure formation may accelerate partic les to 

sufficiently high energ ie s (IColafrancesco fc Blasil 1 19981 : 
lLoeb fc Waxmanl I2OO0I : iRvu et all l2003h . Additionally 
supernovae and galactic winds from cluster galaxies 
can populate galaxy clusters with non- thermal parti- 
cles IjVolk. Aharonian fc Breitschwerdtl 1996). Furthermore, 
powerful AGNs are b elieved to be prominent injectors of 
CRs into the ICM (lEnfilin et all Il997l ; lAharonianl 120021 : 
iPfrommer fc Enfilinll2004l ). 

In this paper we focus on the scenario where a pow- 
erful AGN injects high energy particles into a galaxy 
cluster. Prominent AGNs are often found in galaxy clus- 
ters with short central cooling times of the thermal ICM 
(for a review of these so called cooling flow clusters see 
iFabianl 1 19941 ) and their effect on the ICM can be seen 
in several systems (see iBirzan et al. I 12004 for a sample 
of such galaxy clusters). High resolution X-ray observa- 
tions have revealed bubbles, cavities and weak shocks in 



eral systems (e. 


z. Bohringer et al.lll993l:lBlanton et al.ll2001: 


Schindler et al.l 


2001; iMcNamara et al.l l200ll: iFabian et al] 


20031: Choi et all 20041). Bubbles in the X-ray gas are often 



associated with radio lobes, indicating the presence of rela- 
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tivist i c electrons (lOwen. Eilek fc~ Kassim|[2OO0l ; lFabian et al.l 
20(4 iGitti. Feretti. fc SchindieT 12006t h Very recently, so- 
called cluster-scale AGN outbursts have been found in 
three clusters, all with estimated mechanical energy of 
at least 10 61 erg. Th ese systems are: MS 0735.6+7421 



(|McNamara et al.N2005h. Hercules A (iNulsen et all l2005al ) 
and Hydra A (|Nulsen etaill2005bl ; Twise et al.ll2007h . Large 
scale radio e mission was found in Hydra A with deep VLA 
observations (|Lane et al.ll2004T ). Due to the large amount of 
energy input into these systems by the central AGN, plau- 
sibly in the form of relativistic particles, these systems may 
be promising targets for high (and very high) energy 7-ray 
observations. 

Experimentally, 7-ray astronomy is in a phase of 
rapid development. Several Imaging Atmospheric Cherenkov 
Telescopes (IA CTs) have r ecently been completed: HESS 
Hintoiil 120041). MAG IC (|Lorenj 120041 ) and VERITAS 
Krennrich et alj [20041 ')'). The combination of these VHE 



instruments with the high energy (HE) detector GLAST 
|Thompsonll2004l ). due for launch early next year, will pro- 
vide sensitive coverage of the 100 MeV to 10 TeV energy 
regime for the first time. Here we will argue that these in- 
struments may be close to the sensitivity threshold required 
to detect the high energy electromagnetic signatures of large 
scale AGN outbursts in galaxy clusters. 

In this paper we investigate the high energy luminosity 
of the three known galaxy clusters which host cluster-scale 
AGN outbursts. The relevant properties of these systems are 
given in Table [T] 



2 COSMIC RAYS IN CLUSTER-SCALE 
OUTBURSTS 

The existence of radio synchrotron emission in coincidence 
with X-ray cavities seen in AGN outbursts indicates the ex- 
istence of relativistic electrons within the bubbles and leads 
naturally to the suggestion that relativistic particles may 
support these cavities. This hypothesis must be confronted 
with three key issues: 1. can relativistic particles be confined 
for the required timescales? 2. are the energy loss timescales 
of these particles sufficiently long? 3. is there observational 
evidence for other contributions to the bubble pressure? 

For cluster-scale outbursts the timescales involved have 
been estimated at ~ 10 8 years, considerably longe r than 
for previously identified systems dBirzan et alj 120041 ) . Fig[T] 
shows the relevant energy loss timescales for ultrarelativistic 
electrons in the central cluster environment. We find that for 
the typical B-fields of a few fiG found in the central regions 
of clusters, TeV electrons lose their energy on timescales 
of 10 6 years via inverse Compton and synchrotron cooling 
(see Fig [1} . Only < GeV electrons can exist at the outer 
edge of a bubble of age ~ 10 8 years. The observed radio 
synchrotron emission should therefore provide evidence for 
spectral cooling away from the central AGN. In the case of 
Hydr a A this spectral steepening is clearly seen (|Lane et al.l 
I2004T ). The measurements are broadly consistent with the 
injection of a power-law of electrons (dN/dE oc E~ a ) with 
a ~ 2 which is subsequently cooled above a critical energy 
-Edit ~ 100 MeV to a ~ 3 (resulting in radio synchrotron 
spectral index of ~ 1). By considering continuous injection 
and synchrotron cooling of the electrons over the lifetime of 



the outburst, we estimate that the currently radio emitting 
electrons represent roughly 1/3 of the total energy injected 
in electrons over this period. To estimate the energy car- 
ried by these electrons an estimate of the magnetic field 
strength inside the bubbles is required. Rotation measure 
estimates of the magnetic fie ld outside the radio l obes of 
Hydra A suggest B ~ 30 ^G (jTavlor fc Perlevlll993l ). If the 
B-field within the lobes is similar then the total energy in 
electrons, and in the magnetic field itself, are both close 
to 10 59 erg. As the total energy associated with the Hy- 
dra A outburst is ~ 10 61 erg it seems that at least in this 
system there is missing pressure within the bubbles unless 
the magnetic field strength there is 3> 30^tG. To support 
the Hydra A bubbles solely with magnetic pressure requires 
B ~ 300/iG. Such high fields seem unlikely due to the ob- 
servation of smooth features in the 1.4 GHz emission on 
~ 100 kpc scales. The synchrotron cooling time of the 1.4 
GHz emitting (500 MeV) electrons in a 300 /iG field is such 
that rectilinear propagation at speeds very close to c would 
be required to reach the edge of these feat ures. The missing 
pressure may be ascribed to thermal gas (IGitti et al 1 120071 ) 
or hadronic cosmic rays l|Dunn fc Fabianl [2004) . We note 
that in the past proton dominated jets have been suggested 
l[Celotti fc Fabianl 1 19931 : ISikora fc Madeiskil I2OO0I ) injecting 
hadronic cosmic rays into the ICM during AGN outbursts. 
Cosmic ray protons and nuclei have much less severe en- 
ergy losses than electrons. For a proton to electron ratio of 
only 1/30 these particles could provide the energy required 
to support the bubbles. Furthermore, the escape time of 
< 10 15 eV protons from the central 200 kpc of these clus- 
ters is likely to be longer than the lifet ime of the bubbles 
l|Volk. Aharonian fc Breitschwerdtl Il996l ) . A large fraction 
of the injected CRs may therefore be confined within the 
observed bubbles for the required 10 s years (see Fig. [1} 
but the particles with the highest energies are likely able 
to penetrate the thermal gas which surrounds the bubbles 
(see Fig. [2} • We further note that in cluster scale AGN out- 
bursts the central black hole must accrete a significant frac- 
tion of its own mass in a single activity p eriod and convert 
it very efficiently into me chanical power (|McNamara et all 
120051 : INulsen et alj [2005al lbh. Therefore there is not much 
room for additional radiative energy losses of the relativis- 
tic particles driving the observed shock fronts. 

For Hydra-A it therefore seems likely that thermal gas 
or hadronic CRs dominate the bubble pressursQ- For the 
other objects considered here such detailed radio observa- 
tions do not exist. In general magnetic fields and/or low 
energy electrons can be considered viable alternatives as 
the dominant e nerg y content of AGN jets (see for example 
|Pe Youneir2006l and lDunn. Fabian fc Celottill2006t ). 

Assuming that hadronic cosmic rays are primarily re- 
sponsible for the bubble expansion, two additional inputs 
are needed to calculate the 7-ray flux from these objects. 
Firstly, an estimate of the internal energy in CRs at the 
present day and secondly an estimate of the target density 



1 We note that this conclusion rests on the extrapolation of the 
observed radio spectrum in Hydra-A down to < 400 MeV elec- 
tron energies (i.e. radio frequencies below 74 MHz). Electrons of a 
few hundred MeV have long enough lifetimes to support the bub- 
bles, and could have sufficient energy if there is a second spectral 
component at these energies 
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Table 1. Characteristics of powerful AGN outbursts. The flux F' v is calculated under the assumption of 
scenario A (discussed later). 



Object 


z 


PV 


Age 


mean Density n e 


Bubble Diam. 


vF' v (100 GeV) 






(10 61 erg) 


(10 8 years) 


(10 -3 cm" 3 ) 


(kpc/") 


(10~ 13 erg cm -2 s" 1 ) 


MS 0735.6+7421 


0.22 


6 


1.0 


3 


240/70 


0.7 


Hercules A 


0.154 


3 


0.6 


5 


160/55 


1.2 


Hydra A 


0.0538 


0.9/0.41 


1.4 


5 


210/200 


3.0 




10 10" 
Particle Energy (eV) 



Figure 1. Timescales of relevance to AGN outbursts. For protons 
the mean time between inelastic collisions (solid line) and the 
escape time from the system (dashed line, assuming diffusion as 
given by eq. 9 of Volk et al (1996) with B = 3 fJ.G), are shown. 
Note that this is the escape time from a 100 kpc radius bubble, the 
time required to escape the cluster (under the same assumptions) 
is ~ 10 10 years for a 10 15 eV proton. For electrons the energy 
loss timescale (-E/^) due to synchrotron and inverse Compton 
radiation is shown for three different magnetic field strengths. 
The stars show the mean electron energy contributing to 1 GHz 
synchrotron emission for the three different field strengths. An 
ambient density of 5 X 10 — 3 protons cm -3 is assumed. The target 
photon field for inverse Compton scattering is taken solely as the 
CMB. This is then an upper limit on the cooling timescale, since 
infrared photons from the central galaxy may also act as targets 
for the inverse Compton process. The shaded band indicates the 
apparent ages of the outbursts considered here. 



indicates that the density of bubble material is low compared 
to its surroundings. 

It is probably reasonable to assume an approximate 
pressure balance between the material inside the bubble and 
the ambient gas. Therefore, the external pressure is balanced 
by the total pressure exerted by the sum of the partial pres- 
sures of each component in the bubble: thermal pressure 
(Pth), cosmic ray pressure (Pen), magnetic pressure (P mag ) 
and the relativistic electron pressure (P e i ec ). However, the 
radio data indicate that, in the case of Hydra-A, the mag- 
netic and electron partial pressures are likely less significant 
than the thermal and cosmic ray components. Henceforth, 
we will consider the bubble to be a mixture of only ther- 
mal material and cosmic rays and define Xcr = Pcti/Pext, 
where P oxt is the external pressure on the bubble. In this 
case, the total internal energy content within the bubble is, 



Pint = V 



Pgr 



7th 



7ca 



(1) 



where 7 tn = 5/3 is the adiabatic index of the thermal mate- 
rial, and 7cr = 4/3 is the adiabatic index of the relativistic 
cosmic rays. In this description, the density inside the bub- 
ble is governed by Pcr, and the temperature of the material. 
We note that the case of over-pressurised bubbles Xcr may 
be greater than 1. 

Another useful parameter when discussing the proper- 
ties of bubbles, is the density contrast (77) between the ma- 
terial within the bubble pbub and the ambient material p ex t, 
Tj = Pcxt/pbub- It should be noted that the cosmic rays make 
a negligible contribution to pbub- If we consider a jet that in- 
jects CRs and hot gas at a constant rate, then if this inflates 
a bubble, the average density within the bubble is, 



Pbub 



■ 

W) 



(2) 



where rh is the mass injection rate of the jet, and V(t) is 
the time-dependent volume of the cavity. The enthalpy of a 
slowly inflated bubble is, 



in and around the bubbles. Both these numbers require con- 
sideration of the hydrodynamics of the outbursts, which are 
discussed in the next section. 



E = Z/jott = 



-PcxtV(t) 



(3) 



where Lj e t is the constant jet power and 7 is the effective 
adiabatic index, such that 



3 BUBBLE PROPERTIES 

From the X-ray data we can determine, or constrain, many 
properties of the AGN-inflated bubbles. Observations of the 
X-ray surface brightness profile yield the pressure of the gas 
outside the bubble, the bubble volume, and limits to the 
emissivity of the material inside the bubble. This evidence 



7th 



7" 



7th 



-(1-Xcr- 



7CR 



7CR 



-Xci 



(4) 



Therefore, we see that the volume of the bubble is propor- 
tional to the duration of the AGN outburst. Consequently, 
since the mass injection rate is also constant, so must be 
Pbub- Combining equations ((2} and (|3j and substituting 
Pcxt — ^bPcxtPext /fj, gives the density contrast, 
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v 



7-1 fiL 



jet 



(5) 

7 mfe b r cx t 

which depends on the jet parameters, and the temperature 
of the ambient gas. Assuming a constant mass injection rate, 
r\ increases for greater jet power because the density of ma- 
terial within the bubble is lower. r\ is roughly constant (and 
S> 1) whilst the jet is active. After the jet switches off and 
the bubble rises buoyantly, it will expand to maintain pres- 
sure equilibrium with its surroundings, and the density con- 
trast will fall. We can obtain a simple description of this 
by assuming that the bubble behaves adiabatically. Under 
adiabatic conditions, the specific entropy of the bubble will 
remain constant, thus, 



= constant. 



< 6 > 

Pbub 

Since the temperature of the ICM is roughly constant with 
radius, P cxt oc p cx t, and the bubble is in pressure equilibrium 
with its surroundings, we have: pbub oc p* x 7, an d V oc Pcxt^ 7 • 
Thus, the density of the bubble falls less steeply than that 
of the ICM. 

The argument above describes the slowest possible rate 
at which rj can approach unity. In practise, bubbles proba- 
bly behave far from adiabatically. Analysis of recent hydro- 
dynamic simulations using the FLASH code, suggests that 
the Kutta-Zhukovsky force plays an important part in mix- 
ing the bubble with the ICM and dissipa ting the bubble 
entha lpy over a relatively short timescale (|Pavlovski et al.l 
These simulations show that r\ approaches a value 
of 1, from an initial value of 1000, after rising only a few 
bubble radii. The observed X-ray contrast of the bubbles 
can be used to place limits on the density inside the bub- 
bles for a given temperature. The X-ray contrast is approx- 
imately (ni n t/n cx t) 2 (Tmt/T , cx t) 1/ ' 2 . The uncertainties on the 
X-ray measurements and the geometry of the bubbles are 
such that this quant ity must be less than about ~ 1/3 for 
the Hydra A system 4 Wise et al.ll2007h . Fig. [2] shows this X- 
ray contrast as a function of the relative CR pressure Xcr 
for two different assumed temperatures Tint- 

In the context of this work, it is also important to know 
the density of material in the bubble rim. In the early stages, 
the bubble inflation is likely to be highly supersonic, mean- 
ing that the density of the rim material should be roughly 
four times denser than the ambient material. Since the emis- 
sivity is proportional to the square of the gas density, this 
will lead to the appearance of bright rims surrounding bub- 
bles. However, because the cooling time of these rims is 
short compared to that of th e surrounding g as, t he rims are 
expec ted (and observed, see iBlantonl 12004 and iGitti et al.l 
120071 ) to be cold (fcT < 1 keV). Using the cooling func- 
tion given by ISutherland fc Dopital (|l993l ) we find that the 
time taken to cool from 4 keV (the approximate temper- 
ature of the gas 100 kpc from Hydra-A) to C 1 keV is: 



t coo i » 1.4 x 10 s (no/10~ 



3 )- 1 (feT /4keV) 1 ' 5 years. As- 



suming a compression ratio of four, the rims in these sys- 
tems would have a density 2 x 10~ 2 cm~ 3 , and hence cool in 
about half the age of the observed bubbles. 

2 This force acts to expand the bubble radially, but perpendic- 
ularly to the direction of ascent. As the bubble expands an addi- 
tional component of the Kutta-Zhu kovsky force pushes the b ubble 
down, reducing its ascent velocity l lLandau &: Lifschit j|l995l ). 



4 MODEL 

As discussed in the previous sections we assume here that 
hadronic cosmic rays are responsible for expanding the ob- 
served bubbles. To calculate the rate of proton-proton in- 
teractions, and the secondary particle production in these 
interactions, we apply the parame terisations derived by 
iKelner. Aharonian fc Bugavpvl (I2006TI based on t he SIBYLL 
hadronic interaction model dFletcher et al] [T994). Adapting 
IKelner. Aharonian fc Bugavovl ( 20061 ) (eq~69) , the flux of 7- 
rays from proton-proton (and nucleus-nucleus) interactions 
can be expressed as: 



dNj_ 



Kn e cVe(E 1 



• ^ C{Ep )a pp (E p )^ME,/E p ,E p )^ (7) 



Kn e is the effective cross-section weighted number den- 
sity of target (thermal) protons and nuclei (a thin target 
approximation is appropriate as t pp » tbubble)- Assuming 
a primordial abundance of both CRs and target nuclei we 
find k — 1.15. c is the speed of light, V is the bubble vol- 
ume, e(-E 7 ) represents the 7-7 absorption on the (IR) ex- 
tragalactic background light (EBL). We note that fol l owin. 
recent cons traints from the HESS (I Aharonian et all 1200' 
and Spitzer lDole et al.l (|2006l ) instruments, uncertainties on 
the EBL in the relevant wavelength range are now greatly 
reduced. Here we calculated EBL absorption e(E~,) using 
the wavel ength dependant (z = 0) EBL density given in 
Fig. 13 of lDole et all (|2006l ). ignoring evolutionary effects. 
The function C(E P ) reflects the degree of mixing between the 
target nuclei and the CRs. C = 1 corresponds to the case of 
CRs confined in a uniform medium with electron density 
n e . ^(E^/Ep, E p ) is the energy distribution function for 7- 
rays produced in an a verage interaction of a proton of energ y 
E p (see equation 56 of lKelner. Aharonian fc Bugavovll2006l ) . 
The spectral energy density of protons is assumed to be of 
the form: 



dn p r^—& Ep 
— = kE p exp 

tlll/r) -L^max 



(8) 



With k chosen such that V J 1GeV dn p /dE p = Ecr = 
SP eK tVXcR- Diffusiv e shock acc eleration, either relativis- 
tic (see for example iKirk et aL I l2000l ) or non-relativistic, 
predicts a spectral index of injected particles close to 2. 
We therefore set a = 2 for the purpose of our calcula- 
tions. The predicted peak flux (in the GeV range) varies 
only weakly with the assumed energy range. In contrast, 
the predicted TeV emission depends strongly on the maxi- 
mum energy of the injected particles. The value of at least 
Emax = 10 14 eV (as assumed here) is required to produce 
up to 10 TeV 7-rays. Higher values, common to many pre- 
dictions of hadron acceleration in AGN jets (see for example 
iBiermann fc Strittmatterlll987l ) do not significantly increase 
the > 10 TeV emission, due to the absorption of higher en- 
ergy photons on the EBL. 

The spectrum of the injected secondary electrons is cal- 
culated in a similar way to that of 7-rays. The time evolu- 
tion of the secondary electron spectrum is followed in small 
time steps accounting for synchrotron losses and injection 
via p — p interactions. 
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Figure 2. Left: Relative flux of X-rays and 7-rays as a function of the relative CR pressure Xcr- The X-ray flux from inside the bubbles 
(thin falling lines) is expressed relative to that of an equal volume of gas at the temperature and pressure of the external medium. The 
7-ray flux from the system is shown relative to that expected for a density contrast rj = 1. Two cases are shown: T; nt = T ext (solid line) 
and Tj nt = 4T ox t (dashed line). The dotted horizontal line indicates the approximate range of X-ray contrast excluded by observations. 
Right: The energy dependence of CR confinement within the bubbles (dashed line) and the mixing parameter f for two assumptions for 
Xcr and Ti nt , Xcr = 0.5, T lnt = T ex t (dash-dotted curve, indicated by solid circles in the left panel) and Xcr = 1.0 (solid curve, 
indicated by open squares in the left panel. In both panels the system geometry, external density profile and diffusion coefficient are 
specific to the Hydra A system. Low energy particles do not penetrate far into the external medium and hence encounter mostly low 
density gas inside the bubble. At the highest energies CRs begin to escape from the high density parts of the cluster into the more 
tenuous ICM. Note that the CRs with logio(-E/GeV) = 4 are those primarily responsible for the 1 TeV 7-ray emission. 



The estimation of the quantity ((E p ) introduced above 
is key to evaluating the fluxes from these objects. C(-Ep) de- 
pends on the energy dependent transport of CRs and on 
the density profile of the bubbles. An estimation of the den- 
sity profile in turn requires an estimation of the fraction 
of the bubble pressure provided by CRs (Xcr) and hence 
the remaining pressure provided by hot gas. The left panel 
of Fig [2] shows the relative expected X-ray flux (from in- 
side the bubbles) and the total 7-ray flux (from inside and 
outside the bubbles) as a function of Xcr, for two differ- 
ent assumed internal temperatures. X-ray fluxes above the 
dashed horizontal line are effectively excluded by the ob- 
served X-ray contrast of the bubbles. The CR spatial dis- 
tribution is derived by a numerical transport simulation of 
en ergy dependent diffusion assuming the value suggested 
bv lVolk. Aharonian fc BreitschwerdtJ (|l996l ) eq. 9, combined 
with advection. Note that for the purpose of the diffusion 
calculation a magnetic field inside and outside the bubbles 
was assumed to be the same. In general the diffusion inside 
the bubbles is much more uncertain and could have a sig- 
nificant effect on the observed 7-ray flux. The right-hand 
panel of Fig [5] shows (dashed line) the fraction of CRs con- 
fined within the bubble under these assumptions. Integrat- 
ing over this curve, the total CR energy inside the bubbles 
is 74%. The remaining curves on Fig [2] (right) show C(E P ) 
calculated numerically for the two cases marked on Fig [2] 
(left): Xcr = 0.5, Tint = Text (dash-dotted curve, solid cir- 
cles in the left panel of Fig [2} and Xcr =1.0 (solid curve, 
indicated by open squares in the left panel Fig [2}. For these 
calculations we assume and external density profile based on 



X-ray measurements (see Fig. 4 of lNulsen et ai1l2005br) and 
bubble rims with a compression ratio of 4. It can be seen 
that at high energies (> TeV) the results are rather insensi- 
tive to the internal density of the bubbles, as most emission 
arises from the regions just outside the bubbles. 



5 RESULTS & DISCUSSION 

5.1 Expected broad band emission 

Fig. [3]shows the resulting broad-band spectral energy distri- 
butions calculated for Hydra A and MS 0735.6+7421. Fig.g] 
compares the expected 7-ray emission for the three out- 
bursts considered here. The predicted HE and VHE 7-ray 
flux of Hydra A is the highest of the three clusters. The 
larger energy injected in CRs in the other two systems 
does not compensate for their larger distance (see Fig. |4]). 
The last column of table [T] shows the flux expected around 
10 GeV (Fl) under the assumption that the total energy 
in CR hadrons is equal to 3PV, i.e. Xcr = 1. For Hy- 
dra A, four curves are shown, for different assumptions on 
the distribution of target material and the energy in cos- 
mic rays: (&)Xcr. = 1-0,77 = 1 implying very cold mate- 
rial inside the bubbles, (b)XcR = TO with no gas within 
the bubbles, (c)Xcr = 0.5, r\ — 2 and (d) as for (b) but 
with PV = 4.1 x 10 60 erg rather than the 9 x 10 60 ergs 
used for the other curves. The two estimates for PV in 
Hydra-A come from different approaches, a shock model 
jNulsen et all l2005bj) and summ ing the PV cont ributions 
of individual bubble components (|Wise et alj|2007l ) and can 
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Figure 3. Model spectral energy distributions for Hydra A (left) and MS 0735.6+7421 (right). 7-ray emission arising from the decay of 
pions produced in hadronic interactions, is shown with (solid line) and without (dashed line) the effect of EBL absorption. The inverse 
Compton and Synchrotron radiation of secondary electrons and positrons are calculated assuming magnetic field strengths of and 
0.13 fiG for Hydra A and MS 0735.6+7421, respectively. The CMB is assumed to the dominant target photon field for inverse Compton 
scattering. Radio data are taken from Lane et al (2004) and Cohen et al (2005). 



both be considered valid estimates of the outburst energet- 
ics. We find that EBL absorption in the case of Hydra A 
becomes important above 1 TeV. In the most distant system 
MS 0735.6+7421 EBL absorption is already severe above 
about 200 GeV and the expected 7-ray flux drops dramati- 
cally at higher energies. 

Secondary electrons resulting from the p-p collisions 
will lead to the production of electromagnetic signals via in- 
teractions with magnetic and radiation fields. For comput- 
ing the observational signature associated with secondary 
electrons in Hydra A we use a magnetic field of 30 /iG 
l|Tavlor fc Perlevlll993l ) and inverse Compton up-scattering 
by the CMBR only. We note that the magnetic field of rele- 
vance here is that in the regions where most target material 
exists, i.e. just outside the bubbles. We can compare the 
expected synchrotron emiss ion in the radio band with the 
observed level of radio flux l|Lane et al.ll2004l ) since in prin- 
ciple our model could be tested by the expected synchrotron 
signal from secondary electrons. However, the observed ra- 
dio emission seems to be associated with primary electrons; 
exceeding the expectation for secondaries by a factor ~ 50. 
In the 2-10 keV X-ray band the th ermal emission of the Hy- 
dra A lobes, ~ 10~ erg cm" " 2 s _1 jEdge. Stewart fc Fabianl 
Il992i ) exceeds our predicted secondary synchrotron emission 
by 3 orders of magnitude (see the left-hand panel of Fig. [3j) . 
It therefore appears that in the case of Hydra A the emis- 
sion of secondary particles is completely buried by other 
processes and that the GeV- TeV 7-ray range may be the 
only wavelength band in which the existence of energetic 
hadrons can be probed. 

The situation is somewhat different for 



MS 0735.6+7421. The observed radio spectrum is ver y 
steep, particularly in the outer lobes l|Cohen et all 120051 ) . 
This means that although the low frequency emission is 
likely dominated by primary electrons the observed radio 
flux of the lobes at ~ 1 GHz is very low (~ 2 mj) - of the 
same order as the predicted radio emission from secondary 
electrons. We find that under the assumption that the total 
energy in hadronic CRs is equal to 3PV that a B-field of 
> 130 nG would produce a high frequency hardening of the 
radio spectrum which is not observed (see the right-hand 
panel of Fig. [3}. A value of ~130 nG (the upper limit in this 
scenario) is rather low for the central region of a cluster and 
lies 3 orders of magnitude below the valu e of equipartition 
with relativistic particles (estimated by iMcNamara et ail 
12005I . as 100 pG). We note that radio observations of 
this object at >1.4 GHz are highly desirable to probe the 
existence of secondary electrons (and hence CR hadrons) in 
this object. 

5.2 7-ray observability of these objects 

It is clear from the above discussion that galaxy clusters 
which harbour extraordinarily powerful AGNs are at the 
edge of detectability with current and near future 7-ray 
instruments if the observed bubbles are dominated by the 
pressure of relativistic particles (Xcr — 1). Hydra A is the 
closest cluster-scale AGN outburst known and appears to 
be the most promising target of this kind for 7-ray observa- 
tions. Furthermore, its emission is only moderately affected 
by EBL absorption in the VHE 7-ray regime. Our calcu- 
lations show that, depending on the detailed properties of 
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Figure 4. Predicted 7-ray emission for the most powerful known 
AGN outbursts. Model curves are compared to the nominal sen- 
sitivities of the GLAST, HESS and HESS-2 7-ray detectors, for 5 
year and 50 hour observations, respectively. For Hydra A four 
curves are shown: (a) Xqr = 1.0,7) = 1 implying very cold 
material inside the bubble, (b)Xcn = 1.0 with no gas within 
the bubbles, (c) Xcr = 0.5, rj = 2 and (d) as for (b) but with 
PV = 4.1 X 10 60 erg rather than the 9 X 10 60 ergs used for the 
other curves. 



the source, it may be detectable using the currently oper- 
ating HESS instrument and the upcoming GLAST mission. 
With detections from both instruments it might be possible 
to determine the shape of the spectrum of CR protons in 
this system. Additionally, the extended 7-ray emission ex- 
pected from Hydra A could be resolved by instruments such 
as HESS, but is not sufficiently extended to significantly 
degrade the detection sensitivity with respect to the point- 
source case. A 7-ray detection of Hydra-A would provide 
a clear signature of hadronic CR dominance of the bubble 
pressure. 

The situation is less auspicious for the other two exam- 
ples of cluster scale AGN outbursts. Both systems are more 
powerful than Hydra A but their distances make detection of 
a 7-ray signal more difficult. Furthermore, their larger red- 
shifts compared to Hydra A imply rather severe absorption 
by the EBL, making a >1 TeV detection of these objects 
extremely difficult. It is very likely that these systems can 
only be detected with the next generation of 7-ray telescopes 
(both ground based and space borne). 
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